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We exhibit the multicritical phase structure of the loop gas model on a random surface. 
The dense phase is reconsidered, with special attention paid to the topological points g = 
1/p. This phase is complementary to the dilute and higher multicritical phases in the sense 
that dense models contain the same spectrum of bulk operators (found in the continuum by 
Lian and Zuckerman) but a different set of boundary operators. This difference illuminates 
the well-known (p, q) asymmetry of the matrix chain models. Higher multicritical phases 
are constructed, generalizing both Razakov's multicritical models as well as the known 
dilute phase models. They are quite likely related to multicritical polymer theories recently 
considered independently by Saleur and Zamolodchikov. Our results may be of help in 
defining such models on flat honeycomb lattices; an unsolved problem in polymer theory. 
The phase boundaries correspond again to "topological" points with g = p/1 integer, which 
we study in some detail. Two qualitatively different types of critical points are discovered 
for each such g. For the special point g = 2 we demonstrate that the dilute phase 0(— 2) 
model does not correspond to the Parisi-Sourlas model, a result likely to hold as well for 
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the flat case. Instead it is proven that the first multicritical 0(—2) point possesses the 
Parisi-Sourlas super symmetry. 

submitted for publication to Nuclear Physics B 

SPhT/92-025 



3/92 



1. Introduction and Overview 

1.1. Introduction 

The past few years have seen considerable progress towards understanding theories 
of conformal matter coupled to 2D gravity. Although some insights were due to advances 
in continuum Liouville theory, stunning progress was made in understanding and solving 
discrete models of 2D gravity Initially, the number of such models available was small. 
Due to recent interest in the subject, however, we now have a plethora of such models 
available. Some of these models are known to correspond to conformally coupled matter 
of central charges C < 1, while the continuum limit of many others still remains to be 
understood. Classifying the possible critical behaviour of lattice models of 2D gravity is 
potentially an important pursuit. Because of their interpretation as toy models of bosonic 
string theory we might hope to eventually attack more physical string theories with higher 
central charges and built-in supersymmetry. 

The perhaps simplest such theory is the one matrix model and its multicritical points 
discovered by Kazakov [|I|. They were identified in 0, 0as corresponding to the (2m— 1, 2) 
minimal models coupled to gravity. Even though this identification is beyond any doubt 
correct, it is less clear why the lattice curvature dependent Boltzmann weights of these 
models conspire to give precisely these minimal models. In fact, we only understand 
this puzzle for the first three cases m = 1,2 and m = 3 |. General (p,q) minimal 
models require several linearly coupled matrices ||; actually two are sufficient ||. Again 
a better understanding of why the Boltzmann weights of these multimatrix models result 
in particular conformally coupled theories would be quite interesting. 

An alternative way to introduce C < 1 matter fields onto random lattices is the 
loop gas construction [BJ[I7[. An advantage over the multimatrix approach is that the 
construction is manifestly physical; i.e. it is clear (through the Coulomb gas mapping) 
why a specific critical behavior occurs. All (p, q) minimal models coupled to quantum 
gravity can be constructed in this way. So far a non-perturbative definition of all these 
models has not been given; it is however possible to derive diagrammatic rules signaling 
the existence of a "string field theory" |§. 

The simplest statistical model allowing interpretation in terms of the loop gas is the 
0(n) model ||. The 0(n) on a fluctuating surface is equivalent to a special one-matrix 
model || and thus can be defined in a nonperturbative way. In this simplest model all 
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loops (contractible and noncontractible) are taken with the same weight n. The partition 
function of the corresponding loop gas on a 3-coordinated fluctuating lattice reads 



F(X T N) = jY 2 -2// n #loops e -A 1 ; e -T(total length of loops) ^ ^ 

c/> 3 graphs loops 

where A is the cosmological constant coupled to the volume v = # vertices of the lattice, T 
is the temperature of the loop gas and 1/N is the string interaction constant coupled to the 
Euler characteristics 2 — 2H of the graph. The right hand side of ( |1 . 1| ) represents a triple 
series: It is asymptotic in 1/N 2 and of finite radius of convergence in the fugacities gi = e~ x 
and = e ~ X ~ T corresponding to vertices empty and occupied by loops, respectively. 

For N and n integers this series coincides with the perturbative expansion of the 
vacuum energy of a zero-dimensional N x N matrix field theory || 

VM ft exp{-iVTr[y (M ) + - £ *J - — £ ^M}} (1.2) 

The corresponding Feynman diagrams can be interpreted in terms of surfaces populated by 
nonintersecting loops ( fig. |1|) Then n gives the number of species of loops and the matrix 
potential V(M) specifies the measure in the configuration space of the random surface. 
With the choice V(M) = \M 2 + ^giM 3 the perturbative expansion of the matrix integral 
creates the cluster expansion of the 0(n) model on a 3-coordinated lattice. 

The integration over the matrix variables 3"^, (x = 1, 2, ...,n can be performed imme- 
diatelly and the result is the following one-matrix integral || 

PMexp {-NTtV(M) + - ^(Tre M£ ) 2 } (1.3) 

2 Jo £ 

The method of orthogonal polynomials used in the "ordinary" matrix models is not 
applicable here but the saddle point method in the limit iV — > oo works equally well. 
Indeed, the action depends only on the eigenvalues Pi, ...,Pn of the random matrix M 
and the integral (|1.3| ) equals, up to a constant factor, to 

N 

Y[ dPi g-^EL ^) +Ev j i°8\Pi-p*\-%j: i jog\2Pt>-p i -p i \ (L4) 



i=i 



This last integral is defined for any real value of n and it is known from the flat case that 
criticality occurs for — 2 < n < 2. This is exactly the interval in which the saddle point 
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solution exists ||. In this paper we will restrict our study to the saddle point solution; 
this is sufficient for classifying the possible critical regimes of the model. The stability if 
the saddle point solution and the nonperturbative effects due to instantons are considered 
in 0. 

It is very useful to parametrize n in terms of the Coulomb gas coupling constant g as 

n = — 2 cos Tig (1-5) 

The model can be in a dense phase with < g < 1 or, if gi in V(M) is tuned appropriately, 
in a dilute phase with 1 < g < 2. In both phases the central charge is given by C = 
1 — 6( v /(? — -j=) 2 . Which models may be described within the two phases? Let us agree 
that for our entire discussion in the present work p > q. Then in the dense phase g = | 
and all minimalB models may be obtained; while in the dilute phase g = - and only models 
with p < 2q can be reached. In particular, the p = 2m — l.q = 2 models are not present 
for m > 3. It is therefore natural to put them in by hand and investigate the loop gas in 
the presence of a general potential 

m+l ^ 

y ( M ) = E k 9kMk (L6) 

k=2 

In order to analyze the possible critical behavior it is sufficient to study the theory on a 
disc, i.e. we will consider the loop function 

W(P) = (— Tr— ^) (1.7) 



where the average is taken with respect to ( |1 . 2| ) . For later use we remind the reader that 
the asymptotics of W(P) for P — > oo is given by 

W(P)=p + ^ + ^ + ... (1-8) 

1.2. Overview, Conclusions and Open Problems 

Our main result consists in establishing the phase diagram of fig. [|. In all phases 
the relation n = —2 cos rrg remains valid; i.e. the phases may be thought of as an infinite 



1 Note that the 0(n) model is not minimal at generic values of n, even if the associated central 
charge is in the BPZ list. Special weights have to be assigned to noncontractible loops, as soon as 
the genus exceeds zero or operator insertions are present (see [[?]]). In particular eq.(1.2)may not 
be used to give a nonperturbative definition for the minimal models. 
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number of branches in the complex n-plane (the branchpoints being at n = ±2). The 
Kazakov multicritical points are situated by construction at the center (i.e. n = 0) of 
each phase. For a given central charge between — oo and 1 there are always two points; 
one lying in the dense and the other in a higher phase. These two models related by the 
duality g — > - are however not equivalent. This is particularly surprising for the minimal 



models: There is a dense realization at q = - and a multicritical one at q = -. We 
argue in section 2.2. that the difference is due to a different set of boundary operators. 
From this point of view, then, the infamous (p, q) asymmetry of the matrix chain models 
arises because the chain models always favor one set of boundary operators over the other 
(with the exception of the Hermitian one matrix model with g = |). Another sequence of 
particularly interesting points in the dense phase is generated by the "topological" series 
g = K We demonstrate that these models appear in our lattice approach as theories 
of "boundaries" (i.e. loops) without "worldsheet" I. Our understanding is aided by the 
realization that the loop equation of the model may actually be exactly solved for rational 
values of g. The first member (g = |) of this series is the model of topological gravity 
whose continuum limit is known [|12|. The correct continuum limit of the higher members 
of the series remains to be found. In section 3. we further establish the above diagram by 
investigating the "perturbation" of the (2m — 1,2) models by the loop gas. Note that our 
approach smoothly interpolates between the physical construction in the dense and dilute 
phase and the higher Kazakov points. One may therefore hope to eventually obtain a better 
understanding of the latter. Let us note that flat multicritical 0(n) models have recently 
attracted some attention in relation to a search for multicritical theories of polymers |13[] , 
TA\. Polymers are obtained by taking a n — > limit in the 0{n) model. Dense and 



dilute polymers have been related to a, respectively, broken and unbroken twisted N = 2 
supersymmetry ||13|1 , and it is natural to ask whether some insights into 2D gravity and 
noncritical string theory can be derived from this fact. A further very interesting problem 
is to find out what kind of multicritical polymer theories are obtained by taking ann->0 
limit in the Kazakov backgrounds. A comparison of critical exponents (see 3.3.) indicates 
that our multicritical polymers do not correspond to the ones considered in []T3j . In section 
4. the phase boundaries between the various regimes are studied. A discontinuity at these 
points is detected, resulting in two qualitatively different models at each integer g. The 
first boundary, g = 1, possesses C = 1, while the second boundary, g = 2, has C = —2. 



This way of looking at topological theories was first clearly stated in 
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We demonstrate that the model connected to the m = 3 phase possesses the Parisi-Sourlas 
"target space" super symmetry. The dilute 0(— 2) model is different; only a n — ► — 2 
limit in the dilute phase yields again the Parisi-Sourlas model. The further boundaries at 
g = 3, 4, ... all correspond to "topological" models. Their correct continuum description 
remains obscure for the moment. 



2. The Dense Phase, Revisited 

2.1. The Loop Equation and its Solution 

The dense phase corresponds to setting all couplings {g^} in to zero, except for 
go = 1. All loops are then densely packed on the ip 3 random graphs. The loop equation 
reads 

W 2 (P) = $ 2 ^p { _ P p >) [y'( P ') - nW(2P - P')] (2.1) 

It may be derived from the matrix model ( |1.2|) or by combinatorial methods |7]]. After a 
symmetrization with respect to the reflection P — > 2Pq — P the contour integral can be 
performed using the Cauchy theorem and we arrive at the following functional equation ( 
for the gaussian potential V{P) = \P 2 ) 

W(P) 2 + W(2P -P) 2 + nW(P)W(2P -P) = PW(P) + (2P -P)W(2P -P)-2 (2.2) 

The critical point is located to be at = \/2(2 + n). If we transform P = Pq + P and 
W ( p ) = 2ii^t P + cos n 9 (2P - P)] + W(P) ^simplifies to 

W(P) 2 + W(-P) 2 + nW(P)W(-P) = — ^(P 2 " P*) + 1T—P 2 ( 2 - 3 ) 

2 + n 2 — n 

In [15[,[16| the equation was solved in the vicinity of P*. After introducing a cutoff a we 



blow up the vicinity of P* by the following change of variables 

Po ~ P * = a^A, P ~ Po =az, P °- PR =aM (2.4) 
P* P* — Pl P* — Pl 

Here A is the continuum cosmological constant, z the continuum boundary cosmological 
constant, and M is the boundary cosmological constant induced by the fluctuations of the 
worldsheet geometry. 
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One obtains for the singular part a 9 w(z) of W(P) 

w (z) = -\A g [{z + \/z 2 - M 2 ) 9 + (z- Vz 2 - M 2 ) 9 ] 

= A g M 9 cosh(gT); ( 2 - 5 ) 
z = M cosh(r) 



with 

A = M 2g (2.6) 

A g ==^ — (2.7) 
1 — g sin irg 

The last relation follows if we compare ( |2.5|) with the exact solution obtained for Pq = P* 

m 

- 1 . u 1 - 9 + u 9 - 1 
W = -2 A ° u + u-i 

u + u- 1 1 ' ' 



m 



1 -9 

The inverse Laplace image of w(z) gives the amplitude w(£) for a disc with fixed length 
£ of its boundary 

w(z) = / d£e- ze w{£); 



o 



2 3 / 2 aM 9 1 f 2p^M5-V2£-3/2 e -M£ if £ > j / M 

^ (1 — g) sin 7T(j ' ' ' 

Now (|2.5|), ( pT6|) is quite puzzling for several reasons: 

1. (|2.6|) says that quantum area does not scale like the square of quantum length in 
these models. An attempt at an explanation is made in section 2.2. 

2. For the topological models g = - we see that the string susceptibility is an integer: 
7 S = 1 — -. This is surprising, since it means either that the exact lattice solution is 
analytic in the lattice cosmological constant Po around P* or that we have logarithmic 
scaling violations at these points. Below we will argue that the first possibility holds (in 



accordance with |^8| ). 

3. A subtle but important point is that, for g < |, we have to approach the critical 
value of the boundary cosmological constant as P = Pq + az, not P = P* + az. We do 
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not have a deep explanation for this, but note that it constitutes annother instance of an 
"analytical redefinition" (discovered in ]TT|] ) necessary to get sensible results from lattice 
gravity models. 

It is interesting to realize that ( |2.2| ) may actually be exactly solved for rational g. 
This serves as a check for ( |2.5| ), ( ^.6| )and will help in getting some insight into the second 
problem just mentioned. At the critical point Pq = P* the solution to ( |2.3| )has been known 



already for some time |T7|1 101 ■ ^ mav be parametrized using circular functions or rational 
functions (g = |): 

y/2 1 t 2 P~1 + f 1 



W 



1 - g smitg t 2 P + 1 ,^ 1Q x 

4 , t p 



P = \/l + cos nq — 

l-o V y t 2 P + l 
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Here t is defined on (the double cover of) the Riemann sphere. It is apparent from 
( P-lOp that the Riemann surface of W(P) is algebraic. Furthermore, since the surface 
is uniformized by rational functions of t we see that it has genus zero. P is defined on a 
surface with p sheets. On the lowest (physical) sheet there is one cut from Pj, (branchpoint 
of order 2) to Pr = Pq (branchpoint of order p). The higher (unphysical) sheets possess 
also a cut from Pq to 2Pq — Pl, except for the topmost, which again has just one cut. Now 
it is important to understand what happens if we move away from the critical point: The 
degeneracy at P = Pq will be removed and all branchpoints will be of order 2. However, 
we still have a Riemann surface; it has a finite number of sheets and only square root 
singularities. By a standard theorem we therefore conclude that the solution of (|2.3| )is 
algebraic and thus given by 

W p + hp-^WP- 1 + ... + hi(P)W + ho(P) = (2.11) 

Here the hi(P) are polynomials in P, with Po dependent coefficients. These coefficients are 
determined by expanding ( |2.3|) , (p. 11] ) around P = oo and matching the two expansions. It 
is easy to see that the genus of the Riemann surface will not change off the critical point; we 
therefore know the existence of a parametrization of ( |2.11| ) using rational functions of one 
parameter. Unfortunately we did not succeed in finding a simple such parametrization for 
arbirary g. Nevertheless we are now in a position to get some insight into the second prob- 
lem mentioned above. Given the just presented algorithm for generating exact solutions 
to the loop equation we see that the singularities in the lattice cosmological constant Pq 
are necessarily algebraic for rational g. Logarithmic scaling violations are thus impossible 
and we have analytic behavior in Pq at the topological points g = ~. 



2.2. Bulk and Boundary Operators 

The loops in 2D quantum gravity may be considered as generating functions for the 
local operators of the theory []I9] . Since the loops of the dense phase are clearly different 
from the ones in the standard lattice models it will be interesting to investigate their 
spectrum. Recall that Lian and Zuckerman have calculated the bulk spectrum in a version 
of continuum Liouville gravity. They found for the (p, q) minimal model coupled to gravity 
an infinite number of states with Liouville charges a given by (as always p > q) 

a p + q — k , 
- = ^ fc = l,2,... 2.12 

with the important restrictions k ^ 0(mod q) and k ^ 0(mod p) (7 is the charge of the 
identity operator, k = p — q). In [^0[ ,[plfit was shown that it is precisely these states which 
propagate around the torus. Now, in ||it was found that the torus partition function for 



the dense minimal models coincides with the one found in ||21|| , establishing that the dense 
models do have the expected bulk spectrum. As is well known, the standard matrix model 
realizations of the minimal models contains further operators. One way to see them is 
to look at the scaling dimensions obtained from expanding the macroscopic loops in the 
length t. It is found that there are also operators with charges ( |2.12| )but k divisible by 
p. They were identified by Martinec, Moore and Seiberg ]22[ as boundary operators. The 
first such operator, k = p, has a = ^ and is beautifully interpreted as the length (<^> is the 
Liouville field) 

' ' (2.13) 



For the dense phase models, the spectrum of dimensions found in the loop expansions 
was worked out in jn|. It is easily seen that they correspond to operators with — as in 
( |2.12|) but with k divisible by p missing and k divisible by q appearing. The dense models 
are therefore different from the standard models because they possess a different set of 
boundary operators. This has dramatic effects, like the non-trivial Hausdorff dimension of 
the boundary. Indeed, here the first boundary operator has a = | ^ and may be interpreted 
as the quantum length 

£= le^ (2.14) 



The Hausdorff dimension 2g = 2^ may be directly read off from the last equation. Note 
that the choice of a in ( |2 . 14| ) corresponds to taking the "wrong" branch in the KPZ dressing 
of the identity on the boundary. It appears that the Seiberg rule for choosing these branches 
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does not necessarily hold for boundary operators [ Z5HP2 [. It should also be pointed out 
that the dense minimal models possess p — 1 loops as opposed to the dilute minimal models, 
which have q — 1 loops. In the KP description of the minimal models inspired by the matrix 
chains the dense phase is formally recovered by interchanging the operators P and Q (see 
ID, []5|, [J24| ). This remains however a formal excercise since it is not hard to seeH that 
for the matrix chains the order of P is always larger than the order of Q (except for the 
hermitian one matrix model which is in the dense phase.). Our interpretation of dense 
models as being endowed with "standard worldsheets" but "nonstandard boundaries" is 
especially manifest in the example of the subsequent section. 



2.3. g = l 

g = g is "pure gravity" (C = 0) in the dense phase. And indeed, for n = 1 in the 
maximally dense case (only go = 1 ^ 0) we can integrate out the M matrix in ( |1.2| )and 
obtain 

r i „ i 

(2.15) 



Of course the integration is this trivial only for the partition function and gets harder 
(although it may still be done) if the boundary is present: Remember the loop operator 
is jjTre LM , not jjTte L '^ 1 . It shows, however, that the bulk consists of ordinary (p 4 lattice 
gravity while the boundary is different. 

The exact solution ( |2.11| )to the loop equation (|2.3|) reads 



- (i(P 2 - 6) + P*)W - |(P 2 - 6)1 + 2 -P Q P* = 



(2.16) 



Being of third degree, it may be solved for W(P) to yield (A = Pq — 6) 



W(P) 

h 

3 



1 A 3 


1 


— A2 + - 




27 


v/27 


1 A 3 


1 


A2 - 




27 


s/27 



p Ai 



P\/2(P 2 + 3)P2 _ p4 _ 1 A 2 _ 2^ A | 



(2.17) 



We acknowledge a discussion with M.R.Douglas on this matter. 
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This model is the first nontrivial point in the "topological" series g = ^ and has 
n = — 1. Because of the minus signs one easily proves that all closed string diagrams are 
identically zero because of precise cancellations. The free energy is zero even before taking 
the scaling limit! Open string diagrams, however, survive (see fig. 0). Just as g = -|, the 
model is rather a theory of loops than of surfaces. The solution ( EOT] ) of (273) is 



2 + -P 2 )VT + -(l + P 2 )P- — P 



27 







The equation may be solved to give (A 



W(P) = 
+ 



— P 3 - -(P 2 + 1)P + A/-(P n 2 + 1) 2 P 2 - — A 2 P 2 - — P n 2 P 4 - —A 3 
27 3 K ; V 9 V ; 27 27 27 



+ 



27 



1 P 3 ~ l(PZ + l)P ~ JliPo 2 + l) 2 ^ 2 " i^p2 
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27 



1 -R?_p4 _ 4 A 3 



27 



27 



(2.18) 



(2.19) 



It is manifest that the solution is analytic at Pq = P 2 = 2. 



2.5. g 



In order to study the puzzling question of the singularity in the lattice cosmological 
constant we have analyzed yet annother topological point. The model is less trivial than 
g = | because here n = — y/2 and we do not have precise cancellation of the worldsheet- 
loops away from the scaling region. The exact solution for W(P) is (we introduced the 
abbreviations cr± = 1 ± |v2) 



W(P) = -Ja + (P*-4a- 



+ a.P2-^ V2a. (P-B) y' ( P - ■ P, > ) ( P - P L ) 



(2.20) 



which is, in accordance with (|2.11|) a solution to a bicubic equation for W. The positions 
of the cut are given through P R + P L = 4V2a+P - 2B and P R P L = 4a^(P 2 - 4a_) 2 -^ . 
The auxiliary parameter B is determined by 



S 4 - lV2a + P B 3 + A -a\ (P 2 - 4a_)P 2 - ^(P 2 - 4a_ 







(2.21) 
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It is seenll that any singular behavior in the cosmological constant Pq would have to be 
generated from eq.( |2.21| ). A (most conveniently numerical) study shows however that the 
branchpoints of the physical sheet of (|2.21|) are not located at the critical value P 2 = 4a _ 
(i.e.Po ~ 1-082) but are lying off the real axis at ±0.327±0.633i This confirms the general 



discussion given in section 2.1. and agrees with the picture of [I8"fl. Let us further note 
that, since there are no singularities on the real Pq axis, the lattice model does not undergo 
any phase transition in the bulk cosmological constant, i.e. not even for some Pq < P*. 



3. The Multicritical Phases 

3.1. Solution and Scaling of the Loop Equation 

We will now turn on m — 1 couplings {gk} in order to reach the m th multicritical 
phase. The loop equation is easily generalized to 

W(P) 2 + W(2P - P) 2 + nW(P)W(2P - P) = 



= V'(P)W(P) + V'(2P - P)W(2P - P) + J2J2^+iWj-i[P k - J + (2P - P) k ~ J ] 

fc=U=l 

(3.1) 

where V is the derivative of the potential (|1.6|) . A new feature compared to the dense 



phase equation ( |2.2|) is the appearance of the (as yet) unknown functions {Wj} (see eq.( p..8|) ) 
which depend on all matter couplings {gk},Po- As for the dense phase the linear terms 
may be eliminated through the transformation P = Pq + P, W(P) = 2sin \ ^ [V'(P) + 
cos ng V'(2P - P)] + W(P): 

rn 

W(P) 2 + W(-P) 2 + nW(P)W(-P) = S ^P 2 ' (3-2) 

i=0 

Here the {s2i} are both explicitly and implicitly (through the functions {W^}) dependent 
on the matter couplings. In the m th critical phase the loop function should behave for 
P — » as W(P) ~ P 9 with m — l<g<m. The critical point is found to be fixed through 
the conditions^ 

S = S 2 = ■ ■ ■ = S 2 (m-1) = S 2 m ^ (3.3) 



4 That is because the j-leg functions {Wj} are simple rational functions in B, as ( p.20[ ) shows. 
E.g., Wi is given by P + ^[{V2B 5 - (1 + \V2)P B A - 2a + PgB 3 + (2a 4 + (P^ - 4 CT _)P - (f + 
f V2)P$)B 2 - y/2a%(P% - 4*-) 2 B - (| ^2 + 1)4(P 2 - 4a_) 2 P ] 

5 Note that we have m couplings at our disposal in order to fulfill the m conditions. 
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At that point eq. (|3.2| )may be explicitly solved; the solution with the correct asymptotics 
at P = oo reads 



W(P) 



1 



1 



+ 



2 sin ng sj2 
_Pl_ 



n 



-n* P r 



m-g 




m—g- 



(3.4) 



The position of the left branchpoint Pl as well as the critical couplings {g^} are determined 
from the asymptotics. We see that indeed W(P) ~ P 9 . In view of this solution at the 
m-th multicritical point our arguments about the structure of the Riemann surface of W in 
section 2.1. may be repeated; we thus conclude that the exact solution for rational g = | 
and arbitrary matter couplings is given by 



W q + hg-^W"- 1 + ... + h!{P)W + h (P) = 



(3.5) 



Again the hi(P) are polynomials in P whose degree increases with p; they (as well as 
the functions Wj) are determined by matching (|3.5| ) and ( |3.2|) at P = oo. In order to 
infer the scaling limit of W(P) we however do not need this exact solution; instead we 
simply generalize the arguments of [15 , [16]: upon scaling P = az one concludes from 
( [3 .4]) that W(P) ~ a 9 w(z); thus the dominant piece in (|3.2|) is so which has to vanish 
as a 2g . Parametrizing the proper approach to the critical point (|3.3|) by the continuum 
cosmological constant A one expects to find so ~ a 2g A 9 ; this will give as expected w(z) = 
sinn [(^ ^ V z 2 — A) 9 + (z — \J z 2 — A) 9 ]. The Hausdorff dimension takes on the "classical" 
value 2 for all multicritical phases. 



3.2. Boltzmann Weights on Flat and Random lattices 

The dense and dilute loop gas has been constructed as a direct adaptation of flat 
lattice models to the random case. It is therefore natural to ask whether we can go back 
from our formulation of the multicritical loop gas and deduce how a multicritical 0(n) 
model might look like on flat honeycomb lattices. Generically we would expect to obtain 

]T3"[ . In the case of the 



such a model by introducing interactions between the loops 



first multicritical phase (m = 3) a concrete suggestion may be made, using the procedure 
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of 0. In this case the matrix model ( |1.2| ) generating the loop gas may be rewritten with 
the help of an additional matrix \I/ as 



x 



i' 1 l r~ i n i 11 (3 ' 6) 

exp{-iVTr[-M 2 + |M 3 + -* 2 + i Jf *M 2 + - £ K - E ®l M ^ 

u=l ^ M=l 



which is seen to possess the graphical expansion exemplified in fig. |4|. In addition to the 
loops we have dimers placed on the random "honeycomb" lattice. Their fugacity has to 
be chosen negative in order to perturb away from the m = 2 regime. Note that this 
formulation may be applied to flat lattices and, even though it does not constitute an 
exactly solvable model, could be checked by numerical analysis of the transfer matrix. 
According to our phase diagram in fig. |2| we predict central charges between C = —2 
and C = —7. For the higher phases it is less clear how to proceed. One guess would be 
to complement the dimers by increasingly longer strands (2-chains, 3-chains etc.) with 
alternating signs for their Boltzmann weights. An interesting question from the point 
of view of polymer theory is whether the multicritical matter induces in the continuum 
limit an effective interaction between the loops. An alternative possibility is that directly 
introducing such interactions results in yet another class of multicritical polymers. 

3.3. Geometrical critical exponents for the 0(n) vector model 

The order parameters (the magnetic operators) in the 0(n) model have a simple 
discription in terms of the loop gas. The m-th magnetic operator ^l is represented as 
the source of L nonintersecting lines meeting at a point. The correlation function of two 
such operators can be evaluated as the partition function Fl of a network consisting of 
L nonintersecting lines tied at their extremities, moving in the sea of vacuum loops of 



the 0(n) model | 25[| . The dimension of the star operator ^l can be extracted from the 
dimension of the partition function Fl. This last quantity can be calculated immediately. 

First observe that the L nonintersecting lines cut the world sheet into k pieces with the 
topology of a disk. Let £±, ...,£l be the lengths of the L lines which form the watermelon 
network. As usual, we first sum over all configurations of the world sheet populated 
by vacuum loops keeping these lengths fixed. Then the partition function Fl can be 
represented as an integral over £±, ...II of the product of L loop amplitudes with lengths 

£i + £2,£2 + £3,-Jl + £i 
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dLe~ 2P ^ 



w(e 1 + e 2 )w(e 2 + e 3 )...w(e L + e 1 ) 



By the general scaling argument 

F L ~ A 25 ^-T«f 

where 5l is the gravitational dimension of the operator ^> L . On the other hand the 
loop amplitude behaves asymptotically as w (£) ~ I' 1 ' 9 and therefore 

F L ~ M Lg = K Lg ^ 2v) 

It follows that 



The smallest dimension Ai is positive only in the interval —1/2 < g < 2. Outside 
this interval the propagator of the nonintersecting random walk grows with the distance 
between its two extremities. The fact that two points are connected with a line leads to an 
effective repulsion between them. Such a phenomenon is typical for nonunitary theories. 

4. The Phase Boundaries 

4-1. Exact Solution and Scaling of the Loop Equation at n = ±2 

We will now investigate in some detail the loop gas with n = ±2 in the presence of 
multicritical matter. This will serve as an important check of our ideas. As we discussed 
before we expect to be able to tune to points with g — p integer, which form the crossover 
points from the p-phase to the (p + 1) -phase. Consider again the scaled loop function 




The corresponding flat conformal dimensions are 




w(z) 



-1 



y/ z 2 _ M 2)fl] 



(4.1) 



smug 
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Were it not for the "wavefunction renormalization" factor -A — we would get a trivial (i.e. 

suing ° v 

purely analytic) result for w(z) for g G N. Including this factor and subtracting terms 
analytic in z we instead obtain 



w(z) = ^ [(z + Vz 2 - M 2 ) 9 -{z- a/z 2 — M 2 ) 9 } log \, (4.2) 

7T M 

Upon Laplace-transforming this result, one obtains correctly w(£) = jM 9 K g {Mt\ as one 
expects from Liouville theory JTIJ . We will now demonstrate how to derive this result from 
our lattice model. The algebraic approach clearly breaks down for integer g: the Riemann 
surface associated with W(P) becomes ininitely foliated. Fortunately it is however possible 
to directly solve the saddlepoint equation of the 0(n) model's matrix model formulation 
|[26|| . This equation for the eigenvalue density of the matrix model reads [17 

where V(x) is the potential of the matrix model. It is convenient to change variables to 
A' = \(P Q - x), p' = \(P Q - y), b = |P , p(x) -> \p{\'), gk -> {-\) k ~ 2 9k and rewrite 
as 

r b' 2 , m 

/ PMyrh* =2 E»+i(A'-M fc (n = +2) 

Jo' A -/X 

(4.4) 

/ PW)T!i Jt = 2j29k + i(\'-b ) k (n = -2) 

A further transformation A' = V^. + SA with A = \{a' 2 + b' 2 ) and B = ±(a' 2 - b' 2 ) 
simplifies these equations to 

' fc=i 



f 1 dp^- = 2j]g k+l (VA + BX-bo) k (n = +2) 
r dp ^ ^— = 2 y 9k+1 {JATBX - b ) k (n = -2) 



(4.5) 



y/A + BpX-p ^y/A + BX 
These integral equations may now be immediately inverted [27] to yield 

P (A) = -A v /i_ A 2 f 1 d p-± — 7 =i=f> +1 ( v /ZTs^-& ) fc 

-1 X- Py/l- p 2 f^y/A + BX 



(4.6) 
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for n = +2 and n = —2, respectively. Since p(X) has to be non-negative it is easy to prove 
the positivity condition 



/l . m 

dn-== 9k+i(y/A + Bjl -b ) k = (n = +2) 

/ d » R , ^TT^J29k+i(v / A + B^-b ) k =0 (n = -2) 

J-i V 1 "/" 2 ^ A + B ^t~i 



(4.7) 



Together with the normalization condition for the density J , dfi f p(fi f ) = 1, i.e. 

!/>t™ =1 ( ^ ±2) (48) 

we are thus given two constraints determining the eigenvalue interval [a', b'] as a function of 
the couplings {bo,gk}- Equations (|4.6|) ,(f4.7| ) ,(|4.8|) constitute the exact solution of the n = 
±2 loop equation in the presence of general multicritical matter. We will now investigate 
the singular structure of the solution; specific examples are presented in the next section. 
The critical points are located by setting b' = [IT]]. At the g multicritical point we 
expect the density to behave for A' — > like 

p(X') ~ X' 9 (4.9) 

In order to locate these points we rewrite (|4 . 6| ) for n = +2 as 



p(X) = -Vl-A 2 V r 2l (A + BXY 

[ rn - 1 1 



(4.10) 



1-A 2 ^ t 2i+1 (A + BA) z / dfi 

l — n «/— 1 



S V-i A - VT^l? 

and for n = —2 as 



r m — 3 i 
L 2 J 



p(A) = -Vl-A 2 V r 2i+1 (A + BX) i+1 * 
n z — ' 

i=0 



2 f ' 11 1 

-Wl- A 2 * VA + 5A / d/x- ~^=n= ( 4 - n ) 
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Here the {r 2i , ^2i+i } and {r 2 i+i, t 2 i} are known functions of the couplings {bo, gk} and the 
parameters A, B (but not of A !) whose precise form we do not need at present. The 
remaining integrals in (|4.10|) , (|4.11|) may be expressed in terms of complete elliptic integrals 
of the third kind. At the critical point b' = they become elementary. Transforming back 
to A' ( |4.10| ) may then be written 



- m — 2 1 r m — 1 i 

- 2 J a I 2 J 



p(A') = -^V^A^ £ r 2l X'^ - ^ W 2 ' +1 log ; (4.12) 

In view of fl4.9| ), the possible critical behavior for n = +2 is now transparent from ( f4.12|) . 
We note that it is immediately obvious that we can reach only points with g odd, as 
expected. There are two kinds of critical points: 

1. m = g (we approach the g th point from below.) It is obtained through the conditions 

r = r 2 = ... = r s _ 3 = 

(4.13) 

t X = t 3 = . . . = t g _ 2 =0 tg ^ 

There are g — 1 couplings {gk} and g — 1 constraints, so we indeed have a critical point. 

2. m = g + 1 (we approach the g th point from above.) The conditions are 

ro = r 2 = . . . = r g - 3 = r 3 _i ^ 

(4.14) 

t\ — ts — • • • — tg- 2 — tg — 

Here one has g couplings {gk} and g constraints; ( [4 .14|) defines again a critical point. 
However, one may imagine relaxing the last condition to t g ^ and still be consistent with 
the scaling law ( |4.9[ ), see below. 

An analogous analysis shows that for n = —2 the points with even g are generated. 
At b' = fliTTD gives 



r m — 2 i 
- 2 J 



P(A') 



7ra 



' E J r 2l+ iA /2l+2 -^^ J ^ +2 A' 2 ' +2 log- Y 



/2 



i=0 



Z=0 



a' + Va' 2 - A' 2 



(4.15) 



The first integral in (|4.11| )is actually divergent for b' — > 0. This simply gives the 
condition £q = 0- Again one finds two distinct kinds of critical points: 

1. m = g (we approach the g th point from below.) It is obtained through the conditions 



ri = r 3 
to = t 2 - 



; r 9 -s = 

tg-2 =0 t ff ^ 



(4.16) 
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2. m = g + 1 (we approach the g th point from above.) The conditions are 

n = r 3 = • • • = r 5 _ 3 = r 3 _i 7^ 

(4.17) 

£o — t2 — ■ ■ ■ — tg-2 — t g — 

Here also one has as many constraints as couplings, but it is consistent to relax t g = in 
2.. 

Having found the critical points we will now scale the density in the vicinity of these 
points: Introduce a cutoff a and define ^7 = a( in addition to the usual scaled separation 
of the cuts k' = K — aM. First concentrate on the points of type 2. Approaching them 
along a line in coupling constant space where the conditions ( |4.14|) , ( |4.17| ) remain satisfied 
we obtain from Q4.10|) , ( |4.11|) for the singular limit of p(X') ~ a 9 p(() 



P(0 = C 9 ~ VC 2 - M 2 (4.18) 



a 9 w(z) being defined as the singular part of W{P) = W(P') = —\ dp! p,_ (here 
P' = I(p - p), EL = -a Z ) we derive from (|4~T8|) 



. y _|_ J 7 1 — M 2 

w ( z ) = (-ly-h^Wz 2 - M 2 log M (4-19) 

Some divergent contributions analytic in z, M were discarded in the course of the deriva- 
tion; the general meaning of such terms was elucidated in ||11[| . Eq.( ^.19] )is already almost 
( fl.2| ). For g = 1, 2 they in fact exactly coincide. (Overall wavefunction normalizations are 
ignored.) For g > 3 there are additional terms of the form z l M g ~ x ~' l \/ z 2 — M 2 log . . . in 
( |1.2|) . In view of eqs. (|4.10|) , ( ^.ll|) it is clear how to generate these terms. Approaching the 
critical point we have to tune the [gu] such that the {r2i}, {r2i+i} vanish as the appropri- 
ate power of k' . This constitutes, for g integer, the analog of the "analytic redefinitions" 



of [JTJin the half-integer case. Only along special trajectories in coupling constant space 
do we satisfy the Wheeler-de-Witt constraint. 

It remains to analyze the relation between the cosmological constant A and the scaling 
parameter M. This is done by investigating the conditions ( f4.7| ), ( [4.81 ). The result isi 

A = M 2 log — (o = l) 

8 flM W ' (4.20) 

A = M 2 (g>2) 



6 We have not carried out a detailed proof except for g = 1, 2, 3, 4 (see examples) but strongly 
suspect the validity of our claim for all g. 
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The critical points of type 1 are essentially different from the type 2 just discussed. They 
exhibit logarithmic behavior in the eigenvalue density; it is easily seen that upon scaling 
we do not obtain ( [4.2|) and ( |4.20|) does not hold. 



The general theory of this section will now be applied to the lowest (and most inter- 
esting) values of g. 

4.2. g = l 

This is the loop gas model at C = 1. m = 1 is the right boundary of the dense 
phase discussed in section 2; ti = 1. The critical coupling is 6* = 2. If we introduce 
the cosmological constant as 5q = 2 + 32a 2 A the relation to the scaling parameter is 
A = (M log aM) 2 . The density, however, contains a logarithmic piece, as explained above. 

On the other hand, m = 2 constitutes the left boundary of the dilute phase. The 
parameters in this case read t\ = 1 — 2gibo and r = g\B. Imposing the condition ( f4.14|) , 
i.e. gi = the logarithmic piece in the density is killed and one obtains 

p(X) = 4- Via' 2 - A' 2 )(A' 2 - b' 2 ) (4.21) 

7T&0 



Working out the positivity condition (4/7) gives A = 6q, while the normalization constraint 
yields 

-^-[Aa'E(k') - (A 2 - B 2 )±jK(k')} = 1 (4.22) 
ottoq a 

K(k'), E(k') are the standard complete elliptic integrals of the first and second kind, 

respectively. Considering the limit k' — ► of ( f4.22| )one locates the critical coupling to be 

&o = more importantly, setting 6q = 7^75(1 + a 2 A) one finds a 2 A = |/c' 2 log p, hence 

proving the first assertion in (|4.20|) . Note that the Boltzmann weights of the model are all 



positive, as is expected for a unitary theoryB. 
4.3. g = 2 

Here the central charge is C = —2. m = 2 is the right boundary of the dilute phase. 
The parameters are i to = — ^0 + fi'i^o an< ^ ^2 = 9i- Imposing ( [4.16| )gives gi = The 



7 Be aware of our redefinition of the couplings {gk} following eq. ( |4.3| ). 

8 Note that for g\ = we are in the dense phase (m = 1). The "critical limit" thus corresponds 
to bo — y 0, meaning infinite Boltzmann weights for loops. In other words, it does not correspond 
to a critical theory at all; the lattice version of "C = —00". 
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density contains a logarithmic term. The positivity condition gives 2a'E(k') = nbo while 



the normalization condition turns out to be identical to ( 4.22 ). One thus locates the critical 
point to be at 6q = ft and finds a 2 A = |/c' 2 log A. Therefore the Boltzmann weights are 
positive (a necessary consistency check since we are approaching this point from the dilute 
phase); however, there exists a logarithmic scaling violation A = M 2 log -=jt. 

m = 3 corresponds to the left boundary of the first higher multicritical phase (the 
Yang-Lee phase). One has to = — b +fi'i^o ~~ 92^0, ti = 9\ — 3<72^o an d r i = Bg 2 . Imposing 
= ^2 = it follows 

p{\) = _LAV(a' 2 -A' 2 )(A' 2 -&' 2 ) (4.23) 
tt6 

Positivity gives A = 6q and normalization B = 2bo. Thus 



and it is obvious (6* = 2) that A = M 2 . The couplings are gi = ^f^, g% = gp-. This 
translates into positive weights for the f 3 vertices and negative weights for the f 4 vertices, 
as it must be in the Yang-Lee phase. Surprisingly, the model is exactly identical to the 
"D = —2" theory solved several years ago ||28|| . It possesses a Parisi-Sourlas supersymmetry 
and is most elegantly described by a zero dimensional supersymmetric matrix model. The 
supermatrix is £ = M + 9^ + ty + 9 + 99 A where 9,9 are Grassmann variables, are 
Grassmann valued hermitian matrices and A is the auxiliary hermitian matrix. The action 

S = NTrJ d9d9[±(-^Z)(tz) + l -T? - |e 3 ] (4.25) 

Upon performing the Grassmann integrals and integrating out the auxiliary field one ob- 
tains 

S = ATTr[~(M - gM 2 ) 2 - - 2gM)\ (4.26) 

By looking at the diagrams the action S is generating one finds the graphical expansion 

gV 3Uo^_]_^#(dimers)^_2^#(loops) (4.27) 

{<j2 3 graphs} {dimers + loops} 



One may equally well use a Gaussian propagator. 



21 



where v is the total number of vertices and vq the number of vertices not occupied by 
either loops or dimersS. The dimers and loops are totally self-avoiding (see fig. ^) . Now, a 
comparison of our loop gas Boltzmann weights with (£07D immediately reveals the identity 
of the two models (the two lattice cosmological constants being related through g = 4^)- 
It is interesting to note that the loops are generated by the "fermionic fields" while 
the dimers are associated with the auxiliary field A. What we have shown here is that the 
latter degrees of freedom are absolutely necessary in order to preserve the supersymmetry: 
without them, the critical behavior is different. 

4.4. g>3 

We will only briefly comment on the type 2 critical point m = 4, a g = 3 model. The 
conditions that turn off the logarithmic terms are t\ = 1 — 2gi& + 3(^2^0 — ^9^o = an d 
£3 = 92 — 4<73&o = 0. The third parameter is ro = g\ — "igi^Q + 6^360 + 93A but now we have 
to remember that ro = upon approaching criticality gives ( [4.19|) while we produce (|4.2|) by 
tending r to zero appropriately. It is possible if somewhat tedious to prove a 2 A = |/c' 2 
and therefore confirm (|4.20|) . 

We also carefully investigated g = 4; no new elements or surprises are found. 
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10 Because of the zero mode of the discrete Laplacian ( |4.27| )is of course identically zero for closed 
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Figure Captions 



Fig. 1. A configuration of nonintersecting loops 

Fig. 2. The multicritical phase diagram of the loop gas coupled to 2D gravity. C is the 
central charge, g is the Coulomb gas coupling and n the weight of the loops. 

Fig. 3. An example for a diagram of the topological model (1, 3). 

Fig. 4. A typical diagram in the m = 3 multicritical phase. 
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